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ABSTRACT 

Wavelets are new mathematical objects which act as ‘designer trig functions. To 
obtain a wavelet, the original function space of finite energy signals is generalized to a 
phase-space, and the translation operator in the original space has a scale change in the 
new variable adjoined to the translation. Localization properties in the phase-space can 
be improved and unconditional bases are obtained for a broad class of function and distn- 
bution spaces. Operators in phase space are “almost diagonal” instead of the traditional 
condition of being diagonal in the original function space. These wavelets are applied to 
the squeezed states of quantum optics. The scale change required for a quantum wavelet 
is shown, with Prof. G.M. D’Ariano, to be a Yuen squeeze operator acting on an arbitrary 
density operator. 

1. INTRODUCTION 

Wavelets were created in France less than a decade ago 1 5 when J, Morlet gener- 
alized the phase-space of Gabor 6 by adding a scale change to the frequency^ (^wavenumber) 
axis for applications to geophysical exploration. Grossmann 2,5 and Meyer 3,5,7 immediately 
saw the importance of wavelets for mathematical physics and to deep questions in har- 
monic analysis, respectively. There are a number of review articles 4,5, available today, 
each specializing in different aspects of wavelets. 

In terms of this paper, which applies wavelets to the squeezed states of quantum 
optics, two long mathematics papers are the most important. The author is convinced 
that they will also be the most important for physics, applied mathematics, engineering 
and industrial problems. The two key papers are those of Daubechies , and Frazier and 
Jawerth. 14 Daubechies 13 first constructed a large family of orthonormal bases of compactly 
supported wavelets in L 2 (R n ). Frazier and Jawerth 14 gave a thorough, complete treatment 
of sampled wavelets which is valid both in the classical function spaces and in the modern 
distributional spaces. 

The approach to squeezed states and quantum optics 15 22 will be through the coherent 
states. 23-25 The three main approaches to coherent states are those due to Klauder, 
to Perelomov, 28 and to Onofrio. 29 The Klauder construction starts with an arbitrary rep- 
resentation of a Lie group G on a complex separable Hilbert space H and induces a repre- 
sentation of G on itself with H as a closed subspace of L 2 (R d ,dfi) This yields a subrepre- 
sentation of the regular representation in the sense of Mackey. 30 It works equally well with 
states or frames. The approach of Perelomov starts with a “Little vector and requires a 
multiplier to add enough structure to force projective representation to be unitary I here 
is additional subtlety in obtaining an invariant measure on the coset spaces used to reduce 
G in that G/H i can have an invariant measure dfi { whereas G/H 2 may not. Ihus, the 
choice of a “Little group” or “stability subgroup” is a sensitive issue in the Perelomov ap- 
proach. The Onofrio construction yields a holomorphic representation of the Lie group. At 
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least in simple cases it gives a complexification of the real homogeneous space M = G / II . 
of Perelomov. For additional structure and the proofs see the nice new monograph of 
Kaiser. 31 

2. SQUEEZED STATES OF QUANTUM OPTICS 

The coherent states for each complex number a are generated from the unique, trans- 
lationally invariant Fock vacuum | 0 > using a unitary displacement operator Z?(a) which 

is defined below. Let a, a + be the Bose destruction and creation operator which satisfy 
the canonical commutation relations 





[a, a] = [a + ,a + ] = 0 


and define D as a Weyl-Heisenberg operator, 


( 1 ) 


D(a) : = exp(aa + — a*a) 


(2) 


Then 


| a >= D(a ) | 0 > 


(3) 


is the ordinary coherent state. In terms of an additional complex parameter £, the two- 
photon squeezed states | C, a > of Stoler 15 and Yuen 16 can be generated using the squeezing 
operator 5(C) 


5(C) : = exp((a +2 - C*a 2 ) 


(4) 


through the action 


|C,<*>=S«)I<»> 


= S«)D(a) | o > . ( 5) 

The states generated in Eq. (5) will be called amplitude squeezed states. These coherent 
states satisfy the uncertainty principle but squeeze one side, say time or frequency, 

exponentially. Naively , it would seem that higher order squeezed operators S^ k \k > 2, 
can be defined through the definition 


S {k \(): =exp((a tk -(*a k ) , (6) 

but a neat paper by Fisher, Nieto and Sundberg 32 has shown the matrix element divergence 

< 0 I S w ( 0 I 0 oo , (7) 

for all k 2! This can be interpreted as either non-analyticity of the vacuum or as operator 
domain problems. The task of defining k-photon squeezed wavelets will be relegated to 
future works. Here a new quantum or operator- valued wavelet of D’Ariano and the 
author 33 will be presented. Let A be an observable 

A | a >= a | a > (8) 

where the states | a > give a resolution of the identity 

1 = f dfi(a ) | a >< a \ , (9) 


284 


I 



where dp(a ) is the invariant measure. The generating function of moments of the observ- 
able ,4 in a state whose density operator is p is given by 

< e '“ A >:= Tr[e ioA p] . (10) 

The probability distribution function P(p, a) is defined as 

P{p, a) : =Tr[\a><a\ p] 


and is the Fourier transform of the generating function of moments with respect to the 
measure dp(a ) 

( 11 ) 


< e ,Q>1 > 


= j dp(a)e iQa P(p,a ) 


A filtered Fourier transform with window function 7 (a) for Eq. (11) can be defined natu- 
rally as 

„iad\ . _ J dp{0)e xaa ^j{a)P{p,a) . (12) 


A c-number wavelet transform analogous to Eq. (12) is given by 

W(p, V ,e) : = — p J dp(6) X { 2 -j 1 )P(p,a) • 


(13) 


In the next section, additional discussions of wavelets will be given. 


3. WAVELETS 

For simplicity of exposition, let feL\R\) be a real or complex-valued finite energy 
signal and denote its Fourier transform by f(k). In L 2 ,/ is guaranteed to exist and by 
Parseval’s theorem |/| 2 = |/| 2 with proper normalization. Choose conventions s.t. 


and 


/(*) = 

roc 

/ f(x)e 2n,kz dx , 


/ — OC 


rOO 

f( x ) = 

/ f(k)e~ 2 * ikx dk 
J — OC 


(14) 

(15) 


If supp(f) C (—1/2, 1/2) and feL 2 (R ) 


Inverting Eq. (16) yields 


sin(7r(x — k) \ 

(16) 

7 r(x — k) J 



(17) 
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where x A is the characteristic function of the interval A C R 1 . With discretization the 

mn th coefficient of / (m, n integers) with “window function” g(x — nx 0 ), which is one when 
(x — nx o) is positive, is given by 


C r 


,,n-s: 


O 2nimkozo 


g(x — nx 0 )f(s)dx 


(18) 


Observe that the n-index is a spatial translation of units of x 0 and the m-index is a wave- 
number translation in units of &o- The joint appearance of (m,n) indicates that C mn 

lives in Z d ® Z d , a phase-space. The scale change x — » xjT' plays an important role in 
the Calderon complex interpolation approach. A dilation is a translation in x (space) 
with a scale change in Fourier transform variable k. Calderon 35 published his famous 
reproducing formula in 1964. The conditions required are the following: (i) Let ij> and 

< f be radial, smooth L 2 (R d ) functions whose Fourier transforms ^ and ip in L 2 (R 2 ) with 
support in a set A , 


A := supp{H),p>) = |A; | 0 < Ci <| k |< C 2 < ooj . (19) 

(ii) For each | k 0 

OO 

]T £(2- I k 1)^(2" | k I) = 1 . (20) 

U— — OO 

Let feL 2 (R d ) with /( 0) = 0 and then 


/(*) = i* \mr i * i) . 

(21) 

V 

Let 


9u(k) = f(k)4 3(2" | k |) 

(22) 

which implies that rearrangement of Eq. (21) into 


£ c *M e2 ' i2 '‘ 

S 

(23) 

has an obvious parallel to Eqs. (14-18). Define the quantities 


<PuA x ) : = 2-" / V(2- I 'x - r) , 

(24) 

Vv,r(s) : = 2-^(2 ~ v x-r) , 

(25) 

M x ) : = 

(26) 

and use these with Eq. (23) to obtain 


c a (k) = <p k * f(s • 2*) : =< /, ip ktS > . 

(27) 


Now 
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(28) 


/(* ) = < f k »■ > ’ 

* a 

is a continuous wavelet expansion with wavelet coefficients given by the c s (k) s of Eq. (27). 
The mathematical importance of the wavelet expansion over the Fourier method, is that it 
generalizes to many function and distribution spaces where Fourier analysis is inapplicable. 
The potential physical importance of wavelet methods is to make possible new formulations 
and calculations of physical models. For computation or for experimental signal processing 
the discrete wavelet transform of Frazier- Jawerth called the ^-transformation is required. 
The (^-transform of / is 


EE < /,¥>*,* > 0*,* 

k s 


(29) 


and holds in general function and distribution spaces. The requirements on the L function 
which (fi must satisfy are: 


(?) I (f(x)dx = 0 , (30) 

J — oo 

(«) ¥>(*) = VrA x ) - d ~ ~ ) » ( 31 ) 

/ °° rlk 

I ^( fc ) l 2 < 00 ■ ( 32 ) 


The inverse problem of reconstructing /(x) from coefficients 1,2 ,4 > 5 i can then be reduced 
to the matrix problem 


/(*») 



(33) 


Observe the correspondence of 2" with a scale change and r(r/s) with a spatial transla- 
tion. A major improvement of wavelets over Fourier methods is apparent from Eq. (16). 

Whereas, f(k) has great localization with compact support, /(#) has terrible localization 
properties in x since as x — ► oo 


sm(nx) 1 

7T X | X 


(34) 


The FVazier- Jawerth improvement in localization via the (^-transform is easily seen for 
^(Jfc)eC'o(i?), where the integration by parts 

(2 xix)-f(x)=-j°° (^yk)e-'’ il -’dk , (35) 

can be repeated r-times to obtain 

(2«*)VW = (-I ) r ' (36) 
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Hence, as x — ► oo 


(37) 


i/wi =°(ff) 

The size of c r>3 (f,tp) depends largely on / in a neighborhood of the point (0,^v) with a 

spread (2 _r ) and far from this point the coefficient decays as | x |~ s . The simplest way to 
obtain wavelets is to decompose the space of interest V\ into a closed subspace Vo and its 
orthogonal complement Wo according to the direct sum 


Vi = Vo © Wo 


(38) 


or schematically 


Vi 

/ \ 

Vo W 0 

In order to maintain simplicity let V\ denote either L 2 or t 2 . Let p be a given function in 
Vi which satisfies the relation 


^ h(n)2 1 ^ 2 ¥?(2x — n) , (39) 

n 

where the set of coefficients {h(n)}^ =1 are a collection of constants, called “masking coef- 
ficients” and the 2 1 / 2 factor in front of 9? is for L 2 normalization. If V\ is the closed, linear 
span of all functions {2 1 / 2 * y>(2x — n)} , 


V\ ~ V n {2 1 / 2 v?(2 • -n} 

and V\ ~ 1 2 and is a (useful) special case of Frazier-Jawerth. 
Proposition: If the masking coefficients satisfy the condition 


sup 

n 



— 2 n) 



< A 


v 


for 0 < A < 00 and AeR 1 then the space Vo of Eq. (39) is given by 

Vo = V n p(- - n) 

with Vo ~ £2 and Vo C Vj . 

Proof: Any f(x)eV 1 can be expanded in <p(- — n)’s 

f ( x ) = 'Y^c n p(x - n ) 

n 

= 2 1 / 2 c n h(m)tp(2(x — n) — m) 

n m 

= 2 1 / 2 c n h(m)p( 2x — m — 2 n) 

n m 

Let 


(40) 


(41) 


(42) 


(43) 
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(44) 


where 


and take 


k = m + 2n 

/(*) = 2 1 ' 2 £ 6 *^( 2 * -*) 
k 

bfc — ^ ^ c n h{k 2 n) , 

n 



(45) 


(46) 


The requirement in Eq. (42) suffices for 

M 2 < I 2 ) ’ (47) 

k ^ n ' 

for every ipeV In the event that {h(n)} n=1 e £ 2 , eq. (42) is automatically satisfied and V 0 
is a closed subspace of V\ . 

Question: Does a function if>eV 0 x exist s.t. 


(i) W 0 = V 0 X = Vx {>{ • - n) and 

(ii) {t/>(- — is an orthonormal basis. 


Answer: Yes; Daubechies 13 in L 2 , by Frazier and Jawerth 14 in t 2 •, Besov spaces Sobolev 
spaces, bounded mean oscillation (BMO) spaces and Triebel-Lizorkm spaces. Such a func- 
tion rp is a wavelet and 

{2 I/ / 2 r/>(2‘ / • -n) : v,neZ } , (48) 

is an orthonormal basis for L 2 (R). This reduces the problem to that of finding a finite set 
of masking coefficients {h(n)}^ =1 . The easiest method for finding these coefficients is due 
to Daubechies in Ref. (13). Assume that the set of non-zero masking coefficients is a finite 


set and let <peL 2 s.t. 


m = 1 


(49) 


Given the wavelet expansion 

<^>{x) = ^M n ) 2l/ V(2* -n) , 

take its Fourier transform to obtain 

#fc) = Y J K n Y 2 l " nk/ 2 ^ k l 2 )l ‘ 2}/2 

=: m(k/2)<p(k/2) . (50) 


It is now necessary to show that 


<?( k ) = Jim 


~(n 


m(k/2 : ’) • <f(k/ 2 


">} 


= Y[™(k/2 } ) 

]=\ 


(51) 
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Phi2(x) vs. x 


Phi4(x) vs. x 



Fig; 1 

The second, fourth, sixth and seventh iteration of the Box function for ^(t)- 


makes sense in L 2 (or V\). This suggests a method of finite approximation providing the 
masking coefficients are known. Let r)o(x)eL 2 (Vi ) s.t. 

fjo(Q) = / rf 0 (x)dx = 1 , 

J — oo 

and iterate 

= ^M n ) 2l/ Vr-i(2z - n) , 

n 

to generate p which is a wavelet, but is not the wavelet ^ of Eq. (25), but rather is that of 
p in Eq. (24) instead. In Daubechies nomenclature p is called & father wavelet and if x is 
identified as a “time” variable the dilations (= scale changes and translations) of p span Vo 
which acts as the high frequency, k — w, content of the full space V\. The function p can 
be thought of as a “pixel shape” in Vo as pointed out by Kaiser. 31 Similarly, is called the 
mother wavelet and the dilations of t f? span IV o, which contains the low frequency content 
of V\ . In Figs. 1, and 2 a mother and father wavelet generated by the choice r)o(x ) = b(x), 
the Box function 


(52) 

(53) 


b(x) = 


1, 1/2 < x < 1/2 

0, otherwise 


(54) 







Phi6(x) and Psi6(x) 


Parent Wavelets 



The father wavelet p(x) and the mother wavelet ip(x) for r) 0 equal to the Box function. 


Since xpeWo C V\ , 

%l>(x) = 

ip(k) = m(k/ 2) ■ tp(k) (55) 


where {h(n)}^_j is a set of masking coefficients for Wq and m(-) is a function analogous 
to m(-) in Eq. (47). Let r = n be a translation parameter and observe that finding ip and 
ij> is equivalent to finding two trigonometric polynomials s.t. the 2x2 matrix 


( m(k) m(k) \ 

m(k + r) m(k + r) J 


(56) 


is unitary. It is useful to consider m(k ) as a phase function which partitions by translations, 

\m(k )\ 2 + \m(k + t)\ 2 = 1 . (57) 

To solve for |m(Jt)| 2 = P(| sim| 2 , | cosx| 2 ) treat P as the probability of a binomial process 
with possible outcomes 

pi = | sin .t| 2 

p 2 = | cosa:| 2 . (58) 
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Then m(k) is the square root of P. There are many solutions since P is an even, positive 
polynomial but only one is needed. Then using m, it is straightforward to find m. This 
completes the discussion of simple wavelets. 


For mathematical physics, operators and their expectations are the objects of interest. 
The spectra of operators give infinite dimensional “diagonalizations” in terms of generalized 
eigenfunctions. In a wavelet basis an operator is “almost diagonal” in a sense discussed 
next. 


of T, 


Let T be an operator, / a function in a normed space which is in the dense domain 


s.t. 


and 


Let 


T : f —* Tf 

(Tf)(x)= f K(x,y)f(y)dy 


(0 V<( x ,y)\< 


00 


0 - y I 


dK(x,y) 

1 

dK(x,y) 

c 

dx 

“T 

dy 

_ |x — y| 2 


f( x ) = C )^jk( x ) 

j, k 

where the i/>jjt(x)’s are a wavelet basis. Then the kernel of T can be written as 

K(x,y ) = '^e jk xl> j k(x)rp j k(y) 


(59) 

(60a) 

(606) 

(61) 

(62) 


where the xjjj ^ s satisfy the estimate 

I'M*)! < aTTb 1 + 1 2 ' •* - *ir'“ • («) 

In order to prove that K satisfies the conditions (i) and (ii) split the sum according to 

E = E + E < 64 ' 

3, k j>jo,k 

with jo chosen so that 

2 J ’° < |x - y| < 2 jo+1 . (65) 

Using the decomposition of Eq. (64), the estimate of Eq. (63) and geometrical sums, it 
follows that 

|A"(x,y)| < — T • (66) 

f - y I 

The same techniques yield 


292 


I 



c 


( 67 ) 


and 



( 68 ) 


There are several consequences of conditions (i), (ii) and Eq. (63). One is that such 
operators map L p — » L p for all 1 < p < oo, solving deep, old problems. Another is 
that {^jk} are an unconditional basis for L p , 1 < p < oo. The proofs work because of 
phase-space localization; if two frequencies are well separated their wavelet coefficients 
are small and if two times are well separated their wavelet coefficients are small. The 
localization structure in Eq. (64) is the reason that the disadvantages of “almost diagonal” 
are outweighed by the advantages. 

4. WAVELETS FOR SQUEEZED STATES 

It is clear that in order to define a wavelet for the squeezed states of quantum optics, 
it is necessary to define an operator which changes the scale. This has been accomplished 
in Ref. (33) in a project with G.M. D’Ariano which was initiated at this workshop. 


Let x(') be an analytic function of the observable A defined s.t. 

X(i)|« >= a\a > 


and 


X(A) 


: = J dfi(a)x\a >< a| 


(69) 

(70) 


This function satisfies the relation 

< X(a) > = Tr[x{a)p) 

= j dp(a)x(a)P(p, a) . (71) 

A dilation operator in the Heisenberg picture is defined on an observable A as 

V nf (A) := Ap . (72) 

The time picture is suppressed since no other picture will be used here although Schrodinger 
picture operators are given in Ref. (17). For some observables, A, the dilation operator 
is unitary but D’Ariano 21,34 has shown that there are important operators of quantum 
optics that are completely positive maps, abbreviated CP, and are non-unitary. In the 
unitary case, which is the only case discussed here, all products of operators are preserved 
and 


IV(i) = x{2V(i)} 
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Using Eq. (73), the operator- valued wavelet transform can now be written as 


W{p\ T],e): = 



Tr{V„,(A)p} 



(74) 


Thus, the dilation operator squeezes any state described by a density operator p. In Ref. 
(34) two examples are presented: 

(i) The unitary dilation of one quadrature of the electric field. This case is applicable to 
a phase sensitive amplification. 

(ii) The CP dilation map of the particle number which is applicable to improving noise 
sensitivity in squeezed light signals. 

and analogies of these have been obtained for a quantum wavelet in Ref. (33) with Prof. 
G.M. D’Ariano. 

5. CONCLUSIONS AND OUTLOOK 

The scale change part of the wavelet dilation is accomplished by the Yuen 16 squeeze 
operator. The application of wavelets to quantum optics is an idea with some potential. 
For example, nonlinear modes and mode-coupling using wavelets should prove useful. The 
quantum squeezed wavelet with D’Ariano should be a good candidate for highly dispersive 
biological media. Future work will focus on these ideas. 
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